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Abstract 

The mean free path is an essential characteristic length in disor- 
dered systems. In microscopic calculations, it is usually approximated 
by the classical value of the elastic mean free path. It corresponds to 
the Boltzmann mean free path when only isotropic scattering is con- 
sidered, but it is different for anisotropic scattering. In this paper, we 
work out the corrections to the so called Boltzmann mean free path 
due to multiple scattering effects on finite size scatterers, in the s-wave 
approximation, i.e. when the elastic mean free path is equivalent to 
the Boltzmann mean free path. The main result is the expression for 
the mean free path expanded in powers of the perturbative parameter 
given by the scatterer density. 

1 Introduction 

The mean free path I is an important parameter which allows to char- 
acterize the different regimes in disordered systems H|. When the 
disorder is strong, there appears a localized phase for Kl <C 1, where 
K is the wave number, whereas the diffusive phase is characterized by 
Kl 3> 1. The strength of disorder is thus measured by the parameter 
l/Kl. 

In most theoretical approaches, the randomly distributed impuri- 
ties are often approximated by point-like potentials [||, [I], || . Although 
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this type of potential greatly simplifies the analytical calculations, it 
has been found that this approximation can lead to incorrect results, 
in particular as far as causality violation is concerned p]. Scatter- 
ing theory with finite size potentials does not encounter any problem 
with causality |7]]. Therefore, it is important to use finite size poten- 
tials in order to describe multiple scattering effects. Effects due to the 
finiteness of the potential size are generally neglected even when the 
description is in terms of scattering t-matrices which describe finite 
size potentials f§, [|. As we shall show, contributions coming from 
the momentum dependence of the t-matrix elements and from spatial 
correlations can be neglected. Nevertheless, one has to take care of the 
way corrections due to multiple scattering are calculated. It is indeed 
shown in this paper that the range of the potential plays an important 
role in the calculation of corrections to the Boltzmann mean free path 
in the diffusive regime, where the scatter er density is low. Finally, the 
mean free path is expanded in powers of the scatterer density up to 
the first subleading order in 2 and 3 space dimensions. 



2 Formalism 

We consider the propagation of an electron in a disordered system 
which is described by the following Hamiltonian. 



" = -;£a + V, (1) 

where the potential V = v « 1S the sum of the N randomly dis- 

tributed individual potentials. The simplest finite range interaction, 
which one can imagine, is the hard sphere potential defined by 

~ ,n I oo for |r — RJ < a, 
Vi(r )=v r-R, ) = \ 1 " ' (2) 

10 for |r - Rj| > a. 

Like the point-like scatterer given by a delta function, this type of scat- 
terer has no internal structure (the wave function vanishes inside the 
potential), but it has a finite size given by the range of the potential a. 
The electron propagator can be cast in the form 



G=G + G Vvj }G = G + G TG , (3) 




where Go is the free propagator and T is the global scattering ma- 
trix of the system. This T-matrix is given in terms of the individual 
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scattering t-matrices by the Watson series |T(J 



N 



T =Y^ U + ti^oti + Yl UG tjG t k + ■■■ (4) 

Each term of this operator series represents a sequence of multiple 
scattering. The insertion of the closure relation J |k) (k| = 1 be- 
tween the operators introduces off-shell t-matrix elements. The hard 
sphere off-shell t-matrix elements, (k'|t^|k) with k 7^ k! 7^ y/2mE, 



can be derived from the off-shell barrier potential of finite height [11] 
by extrapolation of its height to infinity. Their expression is given 
in appendix A. The analytic expression, of each term of the electron 
propagator in Fourier space, G{E, k), given by (||), is obtained by 
inserting the momentum closure relation in the scattering series @. 
The momentum integrations are then performed with the use of ex- 
pressions (A.l) and ( |A.2| ) for the t-matrix elements in 2 and 3 dimen- 



sions and the integrals ( |B.1| ) and ( |B.2| ) of appendix B. Throughout 
the calculation, one has to take care of the fact that the scatterers do 
not overlap. This gives rise to correlations between scatterers which 
have to be accounted for when taking the ensemble average of the 
propagator. 

In contrast to the point-like scatterers case where there is no spatial 
correlation, the ensemble average, noted with a bar, of a quantity Q 
reads 

_ \ r N 

q = z / n dKi n - 2a ) Q({^}i=i,- ,n), (5) 

where Rij = |Rtj| = |Rj — Rj| and 

N 



Z = Jf[dRi Yle(R ij -2a) 

i=l j>i 



is a normalization constant, the same as the spatial term of the parti- 
tion function of a 3D hard sphere or 2D hard disk gas. The Heaviside 
step function accounts for the spatial correlations. 

Using results from the virial expansion, especially the Kirkwood 



superposition approximation [12], the spatial correlation function can 



be written as a product of two-body correlation functions g. 

^He(R lJ -2a) = ^l[ g (R l] ). (6) 



z 1 1 v -J > yr 

j>i j>i 
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Then, the average of Q can be approximated by 

Q = / n dR *n^;M{ R *W.«)> (7) 

J i=l j>i 

where the two-body correlation function g is expanded in powers of 
the scatterer density p 

g{Rij) = go(Rij) + pgi(Rij) + p 2 g 2 {Rij) H , (8) 

5o j 5i > 52 • • • being derived from the virial expansion | |12| . In the case 
of hard spheres, go(Rij) = Q(R%j — 2a) describes the fact that two hard 
spheres at position Rj and Rj cannot overlap, and the other terms g\, 
<72 • • • take into account the average effect coming from the presence 
of other scatterers in the system. 

In the case of point-like scatterers go(Rij) = 1 and gi>\{Rij) = 0. 
Notice that using this approximation, i.e. neglecting the spatial cor- 
relations, for the calculation of (k'|tjGg (i£)t 3 -GQ (i£)ti|k) with the 
expressions (|A.2|) and (A.l) leads to a divergence instead of the cor- 



rect result This divergence is due to the convergence conditions 
of integrals ( |B.1[ ) and ( |B.2| ) . It shows that the use of off-shell t-matrix 
elements is restricted to non-overlapping scatterers. This can be un- 
derstood by observing that the existence of off-shell t-matrix elements 
is related to the finiteness of the size of the scatterer. Therefore it 
makes no sense to use these off-shell t-matrix elements in a point-like 
scatterer approximation. 

In order to calculate the average of the propagator analytically, 
one needs an additional approximation, which neglects the correla- 
tions between non successive scatterers which appear in the scattering 
series (|J). Then 

Q = ^ I dRn / II d-Ru +1 g(R iii+1 )Q{{Ri} i=1 .. n ). (9) 
J i=i 

With the help of this approximation, each sequence of the Watson 
series @ containing only distinct scatterers can be calculated. The 
summation of this class of terms is usually referred to as independent 
scatterer approximation [l4| . It gives the first order contribution in 
density of the self-energy. The approximation (|9|) is justified as long 
as the self-energy £ of the average propagator 

G(E, k) = G (E, k) + G (E, k)£(£, k)G(E, k) 
is expanded up to order p 2 , 

£(£,k) = £ (1) (£,k) + £ (2) (£,k) + --- 
0{p) ©(pa) 
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Indeed, corrections coming from the overlapping of two non-successive 
scatterers are at least of the order p 3 , because they correspond to 
scattering sequences containing at least three scatterers. 

In the s-wave approximation, where Ka <C 1, the first order term 
of the self-energy is simply given by the off-shell t-matrix element 

E«(£,k)»p(k|t(£)|k). (10) 

Notice that this expression is the one obtained in the usual indepen- 
dent scatterer approximation. Its imaginary part is related to the 
individual scattering cross section via the optical theorem and leads 
to the Boltzmann mean free path Iq [H], 



K 



( K In 2 Ka 

2 — m 2D, 



(11) 

in 3D, 



Airpa 2 

where the self-energy is taken on the energy shell, E« = E«(£,k) 
with \k\=K = ^M. 

The second term of the self-energy contains the contribution com- 
ing from the resummation of the infinite series of all multiple scatter- 
ing terms involving two distinct scatterers and a further term coming 
from the excluded volume. It appears that this latter term is negligi- 
ble in comparison with the other terms of order p 2 , because of order 
a d where d is the spatial dimension. The two-body scattering term is 
calculated in appendix C and reads 

2m \ K) \ In 3 K'a In 4 K'a J V ; 

where K' = Ke 1 /2 and 7 ~ 0.5772 is the Euler constant in 2D, and 

53(2) ~ ±—L — p 2 J _ + Ka[ lnETa + 31n2 + ln3 + 7-3 \ 

2m K [2 \ 2 J J 

(13) 

in 3D. 



3 Corrections to the Boltzmann mean 
free path 

The mean free path / is defined as the decreasing rate of the average 
propagator 

G(E,T-x f ) oce"^. 
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In order to compare with the case of point-like scatterers ||], we 
consider the deviation of the mean free path with respect to the Boltz- 
mann mean free path Iq. 

I (fKiSpl SfcEM' 3m£< 2 >\ „. 2> 

rr 1 -{-^-— + ^) +0{p) ' (14) 

where T,^' = ^ £W(i?, k)|, „. This formula is easily derived from 
the definition of the mean free path, 

I = — where K 2 = K 2 — (£« + ^). (15) 
23mK 

For point-like scatterers, the term £W vanishes, there is no off-shell 
contribution. 

In 2 dimensions, for low density p and low energy Ka <C 1, one 
obtains 

I TTp ( 1 3 \ 

+ ^Vh^ + h7i^J- ( 6) 

The excluded volume term appears to be of the same order in Ka as 
the off-shell term SW', i.e. O {(Ka) 2 ). These two terms are also of 
the same order in Ka as the first partial wave term which is the term 



for m = ±1 in expression (AT). They are therefore neglected here. 



The leading correction term in (16) is given by the on-shell t- 
matrix and the subleading correction term comes from the first 
term of the two-body series. These corrections lead to a decrease of 
the mean free path with respect to the Boltzmann value. Neglecting 
the subleading corrections in 1/ In 2 Ka, expression (|l6| ) can be written 
in terms of 1/KIq and reads 

h l + ^w- < 17 > 

tO TT Kl Q 

This expression shows that the corrections are not only written in 
terms of 1/KIq but that there appears a logarithmic correction in- 
volving the size of the scatterer. 



In 3 dimensions, the excluded volume and the off-shell corrections 
are also negligible for the same reasons. Contrary to the two dimen- 
sional case, the leading term is not given by the on-shell t-matrix 
£W. This latter is compensated by the first term of the two-body 
series £( 2 ), given by ( p.l|) . It comes out that the leading term of the 
corrections to the Boltzmann mean free path is given by the second 
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term (C.2) of the two-body series. As in 2 dimensions, the mean free 
path is smaller than the Boltzmann expression. 



H-'(^ + ° (a3) )' < 18 > 

The expression (14) for the deviation of the mean free path al- 
ready exists for point-like scatterers in 3 dimensions M, Instead of 
being of the form pa? /K as obtained in (|lg|), the corrections to the 
Boltzmann mean free path for scalar scatterers are proportional to 
p/K 3 . This shows that the behavior of the scalar point-like scatterers 
is very different from that of finite size scatterers. It is much more 
sensitive to small wave numbers. For finite size potentials, the effect 
of two-body scattering is much less important than it is in the case of 
point-like potentials. Written in terms of the disorder strength 1/KIq, 



expression (18) reads 



I 7T 1 . . 

r'-is <19> 

For the scalar point scatterers at resonance, a similar expression is 
found except that the numerical factor is not 7r/4 but 0.375 (see 
eq. (3.14a) from The effect of the disorder is more important 

for finite size scatterers than for point-like scatterers. 



4 Conclusion 

In this paper we worked out the corrections to the Boltzmann mean 
free path of infinite disordered systems composed of hard disk scatter- 
ers in two dimensions and hard sphere scatterers in three dimensions. 
Calculations with finite size scatterers present several differences com- 
pared with calculations with point-like scatterers. One is led to use 
off-shell t-matrix elements. The excluded volume has to be taken 
into account, yielding spatial correlations when averaging over disor- 
der. These two effects appear to be of the same order in Ka and 
are negligible with respect to a more important effect due to multiple 
scattering. They become important when one is interested in non- 
isotropic scattering, i.e. when more than one partial wave has to be 
taken into account. 

To the knowledge of the author, the result obtained in two dimen- 
sions is new and has not been evaluated even for point-like potentials. 
The corrections to the mean free path cannot be given only in terms of 
1/KIq, indeed they contain a logarithmic prefactor depending on the 
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size of the scatterer. In three dimensions, the comparison of our result 
with the result obtained for point-like potentials shows a discrepancy 
in the numerical factor in front of the disorder strength 1/KIq. The 
decrease of the mean free path is more important for finite size scat- 
terers than for point-like potentials. 
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Appendix A : OfF-shell t-matrix 



Following |TT| , the off-shell scattering t-matrix of a quantum particle 
with energy E = takes the form 



fc2 +°° 



(k'|t(£)|k) 2vra ^ e ^ k ,-e k ) 



2mV 

m=— oo 

k 2 - K 2 

7 /2 _ k 2 {kJm(k'a)J m -i(ka) - k' J m - 1 (k'a)J m {ka)) 

+kJ m {k' f a)J m -i{ka) - KJ m (ka)J m (k' a) m ( 7^ > 

H^(Ka) J 

in 2 dimensions and 

<2 fi2 +°° 

<k'|t(£)|k> = a 2 ^(2/ + l)PKcos(k^k)) (A.2) 

1=0 
( k 2 _ K 2 r 

x i k , 2 _ k 2 kji{k' 'a)ji-i(ka) - k' 'ji-i(k' 'a)ji(ka) 

h\+\(Ka)} 

+ kji(k a)ji-i{ka) - Kji(k a)ji(ka)—7—. > 

h\ +) (Ka)j 

in 3 dimensions. The ji and are the spherical Bessel and Hankel 
functions and the J m and H$ are the ordinary Bessel and Hankel 



functions [19]. In 3 dimensions, the angular part is taken into account 
by the Legendre polynomial Pi. 

These expressions enter the calculation of each term of the scat- 
tering series (Q) in the expression of the electron propagator (^) . 



Appendix B : Useful integrals 



In order to derive the expressions of each term of the scattering se- 
ries (Q) in Fourier space, one has to compute products of propagators 
Go with t-matrices. The following integrals are the basic integrals 
appearing in these products. 

The integral used in the calculation of each scattering sequence of 
the Watson series (|j) in 2 dimensions reads 

r+oc k d ~ 1 dk -A- 

J K 2_ k 2 +ie X\ J ^ ka i) J ™( kR ) 

_ n 

= ~ i l Kd ~ 2 II J m> (Kai)H$ (KR). (B.l) 
i=i 

It is valid under the following conditions 

d + m + Ya=i /■*» = 2p p£Z 
R > E"=i a i 

d + Yll=i \m\ > M- 

In 3 dimensions, one gets 

r+oo k d ~ l dk -A- 

_ n 

= ~ l^ijj^Kai^iKR), (B.2) 

i=l 

if the conditions 

d + 1 + 1 + YZU Vi = 2p 

d - 1 + S?=i M» > 1 
are verified. 

Appendix C : Two-body scattering 

The two-body scattering term consists of the summation of all the 
multiple scattering terms with two distinct scatterers 

(iji) + (ijij) + (ijiji) + (ijijij) H (20) 

where (ij) is a symbolic notation for J2i j^i UGotj. 
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In 2 dimensions, the average of the first term of this series is the 
leading term and can be calculated exactly [15]. 



,(2) 



Siirp 1 



n 2 



h 2 

--16iirp 2 a 2 

2m 



( Mz) \ 


3 


WW 


L 


( Mz) N 


3 




)[ 



r(l)/ 



where z = Ka. When Ka — > 0, expression (12) is recovered. The 
average of the other terms appearing in the series (20|) can be evaluated 

Hq (Kr) I . They lead to 

an infinite series of powers of t — Jo ^ 
energy. 

8in n h 2 



1+^ In feT 



for the self- 



S (2) 



J fC 2/9 2m 



n>2 



(21) 



with a n ss (|)""!- 

In 3 dimensions, the geometrical series is first resummed before 
averaging. Then the order p 2 term of the self-energy is given by the 
average of this resummation and reads 



2m K 



r 2 drh { +\Kr) 2 (l- J^L jo (Kr)h { +\Kr) 
ftp (z) 



2a 



h ( +) (z) V ' 



The integral is decomposed into three parts, noted I^ 2 \ and 
The first contribution represents the scattering sequence (iji). 



/« = 1 

~ z 2 



dx e 



2izx 



'lizx 



2z 3 2z 3 z 2 ' 
The second contribution is the scattering sequence (ijij) 

sin ze~ 



j(2) 



-e 3lzx sin zx 



1 / ITT 

In z + 4 In 2 + 7-1 



(C.l) 



(C.2) 
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In the limit z — > 0, the remainder of the series can be evaluated by 
taking z = 0. This gives 



1 



/(3) = _ 



dx 



z 2 J 2 x(x + 1) 



ilnl. 



Combining these three expressions, one obtains expression (13) for the 
self-energy. 
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